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dance without the usual approximation of partial wave expansion. 
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1. Introduction 



The quest for identifying the nature of the dark matter (DM) in the Universe continues |T|, 
|2|. It is generally believed that most of the DM is made of some hypothetical weakly- 
interacting massive particles (WIMPs). From the particle theory point of view, a commonly 
considered candidate for the WIMP is the lightest neutralino under the assumption that 
it is the lightest supersymmetric particle (LSP). In most approaches the LSP is stable due 
to an additional R-parity ||. The neutralino, being massive, often provides a sizeable 
contribution to the relic density. In addition, the requirement that the neutralino, or 
some other stable particle relic, does not "overdose" the Universe, often provides a strong 
constraints on a supersymmetric model. The two other robust candidates for the LSP and 
cold dark matter (CDM) are the axino (superpartner of the axion) jjj] and the gravitino 0]. 

Continuing improvements in determining the abundance of cold dark matter, and other 
components of the Universe, have now reached an unprecedented precision of a few per 
cent ||. In light of this, one needs to be able to perform an accurate enough computation 
of the WIMP relic abundance, which would allow for a reliable comparison between theory 
and observation. 

The literature on the relic abundance of the neutralino is vast and still growing. (For a 
comprehensive review, see Ref. P].) A brief, and by no means complete, account of major 
developments, can be summarized as follows. The original paper by Goldberg |7[] consid- 
ered the neutralino in the photino limit and pointed out the strong constraints from its 
relic abundance. This was soon followed the first analysis by Ellis, et al. [||] and Krauss ||, 
of the general neutralino case. Several other early papers subsequently appeared with more 
detailed and elaborate analyses. In particular, Griest |l(]] was first to compute in detail the 
annihilation into the ordinary fermion-pair (//) final states through the Z-exchange, and 
later Griest, Kamionkowski and Turner [jllj conducted the first more complete analysis of 
the general neutralino case into WW, ZZ and Higgs-pair final states. The Higgs contribu- 
tion to // was first computed in Refs. [12, 11]. Olive and Srednicki [13] considered all the 
annhilation channels but only in the limit of the pure gaugino and higgsino cases where 
several important resonances and final states are absent. Drees and Nojiri |14] computed 
a first complete set of expressions for the product of the cross section times velocity using 
the helicity amplitude technique. When expanded in the nonrelativistic limit, these give 
expressions for the first two coefficients of the partial wave expansion. 

In the early papers the partial wave expansion of the thermally-averaged product of 
the neutralino pair-annihilation cross section and their relative velocity, (av) ~ a + bx 
was used in most cases. The method is normally expected to give an accurate enough 
approximation (about few per cent) but only far enough from s-channel resonances and 
thresholds for new final states, as was first pointed out by Griest and Seckel [15] and 
further emphasized in Refs. [16, 17, 18, 19|. In particular, it was shown in Ref. [17] that, 
because of the very narrow width of the lightest supersymmetric Higgs h, in the vicinity 
of its s-channel exchange the error can be as large as a few orders of magnitude. Partial 
remedies were suggested, for example in [17, 20], by numerically integrating the full cross 
section near resonances only and by matching this with the expansion-based calculations 
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further away from poles. Such methods are not fully satisfactory since they do not include 
interference terms. A recent detailed analysis showed that in the case of the often wide 
s-channel resonance caused by the exchange of the pseudoscalar A, the expansion produces 
a significant error over the range of neutralino mass which can be as big as several tens of 
GeV. Furthermore subdominant channels and often neglected interference terms can also 
sometimes play a sizeable role. 

A formalism for computing the relic abundance also became refined. In particular, 
the effect of replacing the usually assumed common heat bath for both annihilating par- 
ticles [|]J by a more accurate treatment of involving two separate thermal distributions 
was considered by Gondolo and Gelmini [16| and by Srednicki, et al. [21|. In practical 
terms the numerical difference is usually negligible when one does the usual partial wave 
expansion. The first coefficient a is universal while the second ones b differ by 3/2a, where 
usually a <C b. Gondolo and Gelmini ju]] further derived a very useful compact expres- 
sion for the thermally-averaged product of the neutralino pair-annihilation cross section 
and their relative velocity as a single integral over the cross section, as we will see below. 
Gondolo and collaborators next developed a Fortran code DarkSusy [22| where the relic 
density of neutralinos is numerically computed without using the partial wave expansion 
approximation . 

An additional effect of reducing the relic abundance of WIMPs through co-annihilation 
was first pointed out by Griest and Seckel |jl5|| . In some cases there may exist some 
other states which are not much heavier than the stable WIMP and may therefore be still 
present in the thermal plasma around the WIMP decoupling. In the framework of minimal 
supersymmetry with the lightest neutralino LSP with a significant higgsino component, 
the co-annihilation with the next-to-lightest neutralino and the lightest chargino is often 



important [23, £4|]. Other cases of interest involve neutralino annihilation with the lighter 
superpartner of the r-lepton ]25| and with lighter stop [p6| . 

In this paper, we will present a full set of exact analytic expressions for the cross section 
of the neutralino pair-annihilation in the general Minimal Supersymmetric Standard Model 
(MSSM) for a general neutralino case. From the point of view of low-energy supersymmetry, 
the most natural choice for the LSP and CDM is a nearly pure gaugino (bino) as was first 
shown by Roszkowski Iffifl . Remarkably, just such a case of the [electrically neutral] LSP 
naturally emerges in most case in the Constrained MSSM j28| , [l7| , [29|| . Nevertheless, in 
our analysis we will make no simplifying assumptions about the neutralino, nor will we 
assume the degeneracy of the left- and right-sfermion masses. We will include all tree- 
level final states and all intermediate states. We will also keep finite widths in s-channel 
resonances. We will only neglect possible CP-violating phases in the SUSY sector. We 
will also not consider the effect of co-annihilation here but will address it in a subsequent 
publication. A complete set of expressions presented here does not rely on the partial wave 
expansion, includes all the terms and is valid both near and further away from resonances 
and thresholds for new final states. 

Some of the results presented here are not new but we include them here nevertheless 
in order to provide a complete and self-consistent reference containing the full set of exact 
expressions. In particular, the cross sections for the neutralino pair-annihilation into the 
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SM fermion-pair (//) final states were first computed in Ref [10] and for the WW, ZZ 
and Higgs boson-pair final states in Ref. [11]. 

Given the complexity of analytic expressions presented here and in other papers and 
due to often different conventions used, it is not doable to perform a real comparison with 
the literature. Instead, we have performed a numerical check of our cross section with 
the results obtained by using DarkSusy S|. We have found, for the same values of input 
parameters, an impressive agreement, at the level of a few per cent, for all the annihilation 



channels, which we find reassuring. (Recently another numerical code has been derived [30] 
and also numerically agreed with DarkSusy.) 

While the exact analytic expressions presented here are applicable both near and away 
from special cases where the method of partial wave expansion fails, sometimes one may 
find it less CPU-time consuming to use the latter one. Starting from our exact cross 
sections we have therefore derived a complete set of expressions for the usual first two 
terms of expansion for all the dominant channels. We will present them here as well. 

The annihilation into // is often dominant. However, other final states can also play 
important role, depending on the case. In a previous paper [19| we performed a detailed 



numerical comparison of the relic abundance computed using the exact formulae with the 
one obtained using the expansion formulae, for all the channels, including subdominant 
ones. Our analysis confirmed that the expansion gives highly inaccurate results near res- 
onances and new thresholds. We also showed that very far from such cases the error is 
typically rather small, of the order of a few per cent. However, we found that, because of 
the existence of several resonances (Z and the Higgs bosons), the expansion produces large 
errors, compared to an exact treatment, over a sizeable range of the neutralino mass m x , 
even of a several tens of GeV. In other words, the widely used method of expansion may 
lead to significant errors in a sizeable fraction of the neutralino mass. 

The plan of the paper is as follows. In Sect. |2| we review the formalism for computing 
the relic density that we employ. In Sect. || we introduce the relevant ingredients of the 
MSSM and list all the neutralino pair-annihilation channels. Explicit expressions for the 
annihilation cross secion are given in Sect. |4[ In Sect. |^ we discuss expansion and provide 
a list of formulae for the first two coefficients in the case of equal-mass final states. In 
Sect. ^ we summarize our work. Appendix A contains a complete list of Lagrangian terms 
and couplings which are used in the paper while in Appendix B we provide expressions for 
several auxiliary functions used in the text. 

2. Calculation of the Relic Density 

The relic abundance of some stable species \ is defined as Q x = p x /p cr it, where p x = 
m x n x is the relic's mass density, n x is its number density, p cr u = 3Hq /8ttGn = 1-9 x 
10~ 29 (h 2 ) g/cm 3 is the critical density and Gn is the gravitational constant. (For a review 
of relic density calculations, see, e.g., Refs. Jl|, §].) The time evolution and subsequent 
freeze-out of n x in an expanding Universe are described by the Boltzmann equation 

±j& = - 3H n x - (av M0l ) [n\ - (n^) 2 ] , (2.1) 
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where n x q is the number density that the species would have in thermal equilibrium, H(T) is 
the Hubble expansion rate, cr(xx ~~ ^ an ) denotes the cross section of the species annihilation 



into ordinary particles, VM0I is a so-called M0ller velocity [16] which is the relative velocity 
of the annihilating particles, and (CUM01) represents the thermal average of <7i>M0i which will 
be given below. In the early Universe, the species x were initially in thermal equilibrium, 
n x = n x q - When their typical interaction rate T x became less than the Hubble parameter, 
T x < H, the annihilation process froze out. Since then their number density in a co-moving 
volume has remained basically constant. 

The thermally-averaged product of the neutralino pair-annihilation cross section and 
their relative velocity {(TVM01) is most properly defined in terms of separate thermal baths 
for both annihilating particles [16, ^] 



{aVM0l) {T) = J<P Pl d*p 2 e-ei/T e -E 2 /T ( 2 - 2 ) 

where p\ = {E\, p x ) and P 2 = (-E2, P2) are the 4- momenta of the two colliding particles, and 
T is the temperature of the bath. (Note that one often uses another definition of (<7t>M0i) 
which involves a single thermal bath for both neutralinos. Compare, e.g., Refs. Jl|, [l4|, pfl. 
The numerical difference between the two formulae is usually rather small.) The above 
expression can be reduced to a one-dimensional integral which can be written in a Lorentz- 



invariant form as 16] 



(av M0l )(T) = 1 — dsa^is-imD^K, , (2.3) 

where s = (p\ + P2) 2 is a usual Mandelstam variable and Ki denotes the modified Bessel 
function of order i. In computing the relic abundance one first evaluates eq. ([^) and then 
uses this to solve the Boltzmann eq. (|2.1|). 

There are a number of methods of solving eq. (|2.1|). One often used, approximate, 
although in general quite accurate (for a recent discussion see Ref. |19[|), solution to the 
Boltzmann equation is based on solving iteratively the equation 



I 7 1 = b (5\S<^' > ^)4 /2 ). (") 

where g* represents the effective number of degrees of freedom at freeze-out {^/gl o± 9). 
Typically one finds that the freeze-out point xt = Tt/m x is roughly given by 1/25-1/20. 
One usually introduces J(xf) defined as 

r x f 

J(x f )= / dx(av M 0i)(x), (2.5) 
Jo 

where x = T/m x . 

The relic density at present is given by 



1.66 / T, 



3 



^ = is W T ^iwr <2 - 6) 
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where Mpi = l/y/Gjf denotes the Planck mass, T x and T 7 are the present temperatures 
of the neutralino and the photon, respectively. The suppression factor (T x /T 7 ) 3 ra 1/20 



follows from entropy conservation in a comoving volume [31]. 
3. WIMP Annihilation in the MSSM 

In this Section we introduce the relevant parameters and definitions. We will be working 
in the framework of the general MSSM. (For a review, see, e.g., Ref. §. We follow the 
conventions of Ref. P^j .) The lightest neutralino is a mass eigenstate given by a linear 
combination of the bino B, the neutral wino W® and the two neutral higgsinos and H® 

X=Xi = N U B + N l2 W$ + N 13 H$ + N lA Hl (3.1) 

The neutralino mass matrix is determined by the U(1)y and SU(2)l gaugino mass parame- 
ters Mi and M2, respectively (and we impose the usual GUT relation Mi = | tan 2 0y/M2), 
the Higgs/higgsino mass parameter fj,, the usual weak angle 0w and tan/3 = vt/vb - the 
ratio of the vacuum expectation values of the two neutral Higgs fields. 
The neutralino mass matrix is given by 



M x o = 



( Mi — mz sin 6w cos (3 mz sin 9w sin (3 ^ 

M 2 mz cos 9\y cos j3 —mz cos 9\y sin j3 

— mz s'mOw cos /? mz cos Ow cos f3 — /i 

y mz sin ^ sin f3 —mz cos ^vk sin (3 —fj, / 



(3.2) 



The neutralino mass matrix is diagonalized by a unitary matrix N 



N*M. x aN 1 = diag(m x o, m x o, m x o, m x o). (3-3) 

In the absence of possible CP violating phases, one can choose a basis such that the mixing 
matrix ./V is real, in which case some of the neutralino masses will in general be negative. 
The chargino mass matrix is given by 

M x± =( r M2 yWsin/A 
\ v 2mw cos (3 fi I 

The chargino mass matrix is diagonalized by two unitary matrices U and V 

WM^V- 1 = diag(m v± ,m v± ). (3.5) 

There are two neutral scalar Higgs bosons h and H , a pseudoscalar A plus a pair of 
charged Higgs . (We will typically suppress the Higgs charge assignment except where 
this may lead to ambiguities.) 

Other relevant parameters which determine the masses of scalars and various couplings 
are the squark soft mass parameters mg, mjj and mo, the slepton soft mass parameters m^ 
and mE, and the pseudoscalar mass m^. We also include the trilinear terms A, (i = t, b, r) 
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of the third generation which are important in determining the masses and couplings for 
the stop t\ t 2, sbottom and stau states, respectively. 

In general the flavor-violating sfermion mass-squared (6 x 6) matrices are given by 



, 2 _ ( m| L (/) m 2 LR (f) 
m 2 RL (f) m 2 RR (f) 



M 2 7 = [ ^ L ) J J I > (3-6) 



with 



m 2 LL (f) = M~ L + MjMf + m| cos 2(3 (T 3 j — Qj sin 2 9w) 1, (3.7) 
m RR (f) = Mj+M f Mj + m 2 z cos2(3Q f sin 2 6 W 1, (3.8) 



+ M f^f 

Afj (a^ - ficot f3\ for T 3/ = +1/2 



<n(f) = < I > ( , (3-9) 

1 Mj (Aj - /itan/3j for T 3f = -1/2 



= m 2 LR (f), (3.10) 

where / and / denote here different types of (s)fermions (up- and down-type (s)quarks, 
charged and neutral (s)leptons). For (s)neutrinos only the LL part of eq. (|3.6[) should 



be taken. M~ and M~ denote 3x3 soft SUSY-breaking sfermion mass matrices, and 
jl fR to ' 

Mf denotes here a 3 x 3 fermion mass matrix. Finally, Aj is a scalar trilinear coupling 
matrix of the same dimension while Qf and T^f are the respective electric and isospin 
charges. All the interaction terms and couplings that we will need below are summarized 
in Appendix A. 

In the MSSM, the neutralino LSP's can pair-annihilate into a number of final states, 
if kinematically allowed. A complete list of all tree-level two-body final states is given in 
Table [j]. We only neglect two-body loop processes into final state photon pairs and gluon 
pairs because they are always subdominant in computing the relic density Q. We also 
neglect three-body final states since they are unlikely to be competitive with two-body 
final states. They were shown to dominate in the higgsino case just below the WW and 
ti final states p3| ] but in such regions neutralino co-annihilation with the lightest chargino 
and next-to-lightest neutralino reduce the relic density to very small values anyway. 

The channels WW, ZZ, it, W^H^, Zh, ZH, Ah and AH, are not s-wave suppressed, 
and, once kinematically allowed, can give dominant contributions. But even the s-wave 
suppressed channels // (/ ^ t), hh, Hh, HH, AA, H + H~ and ZA can play some role, 
especially if the other channels are not yet kinematically allowed. This in particular is the 
case with the light fermion-pair final states for which the cross section is suppressed by the 
square of the corresponding fermion mass but which are always kinematically allowed and 
often dominant. 



4. Exact Expressions 

We now proceed to present a full set of exact, analytic expressions for the total cross section 
a (xX ~* all) for the neutralino pair-annihilation processes into all allowed (tree-level) two- 
body final states in the general MSSM. We have included all contributing diagrams as well 
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Process 


Exchanged particles 


s— channel 


t— and u— channel 


XX ^ /m 


h ff 


v° 

Xi 


vv v ffff 

XX ^ nn 


A ff 


v° 

Xi 


XX ^ nri 


h ff 


v° 

Xi 


vv v A A 

XX ^ AA 


h ff 


v° 

Xj 


XX ^ < lJi 


4 7 


Xi 


vv \ H A 

XX —* nji 


fi ff 


Xi 


w — v ff + ff- 


h H 7, 


Xfc 


vv — v ff T 
A A ^ vv n 


ft /7 4 


Xfc 


AA ? Zy "' 


A,Z 


X? 


xx ^-ff 


A,Z 


X? 


xx -> #4 


h,H 


X? 


xx - 


h,H,Z 


Xfc 


xx -> 


h,H 


X? 


xx -> // 


h,H,A,Z 


fa 



Table 1: A complete set of neutralino pair-annihilation channels into tree-level two-body final 
states in the MSSM. The indices i, k, a run as follows: i = 1, . . . , 4, k = 1, 2 and a = 1, . . . , 6. 



as all interference terms and kept finite widths of all s -channel resonances. We have made 
no simplifying assumptions about sfermion masses although we assumed that there are no 
CP violating phases in SUSY parameters. 

To start with, it is convenient to introduce a Lorentz-invariant function w(s) plj 



w(s) = ~ J dLIPS|.A(xx-all)| : 



(4.1) 



where |*4(xx ~~ *■ an )| 2 denotes the absolute square of the reduced matrix element for the 
annihilation of two x particles, averaged over initial spins and summed over final spins. 



The function w(s) is related to the annihilation cross section <r(s) in eq. (|2.3| ) via [ 18 1 



w(s) 



s(s — 4m?,) <r(s). 



(4.2) 



Since w(s) receives contributions from all the kinematically allowed annihilation pro- 
cess XX ~^ /1/21 it can be written as 



W ^ = 32T^ cd ( s ~ ( m h + m h?) Pf(s,m fl ,m f2 )w flf2 (s) 



hh 



(4.3) 



where the summation extends over all possible two-body final states /1/2, mf x and rrif 2 
denote their respective masses, and 



Cf if /l(2) = /(/) 

1 otherwise, 



(4.4) 
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where cj is the color factor of SM fermions (at = 3 for quarks and cj = 1 for leptons). The 
kinematic factor (3f is defined as 



P f (s,m h ,m h 



1/2 



m 



1/2 



(4.5) 



In the CM frame, which we choose for convenience, the function Wf 1 f 2 (s) can be expressed 



as 



1 f +1 

w /i/ 2 ( s ) = 2 J 1 dcos °CM \A(xX /i/2)| 2 , 



(4.6) 



where 9cm denotes the scattering angle in the CM frame. In other words, we write |-4.(xx — * 
/i/2)| 2 & s a function of s and cos 6cm, which greatly simplifies the computation. 

We will follow Table [l] in presenting explicit expressions for Wf 1 f 2 (s) for all the two-body 
final states. All the couplings are defined in Appendix A. All other auxiliary functions, are 
listed in Appendix B. Feynman diagrams corresponding to all the annihilation channels 
are given in Chapter 6 of Ref . M . 



1- XX — ► hH 



This process involves the s-channel CP-even Higgs boson (h and H) exchange and the t- 



and u-channel neutralino (Xi 



w h H 



- 1 , . . . , 4) exchange 
w hH + w hH + w hH 



(4.7) 



CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 
hH 



(jhHr (jXXr 



y 

s — m 2 + iT \ m r 

r=h,H 



4 m, 



(4.8) 



neutralino (%?) exchange: 



w 



(x°) 



hH 



m v o mnllf + m x mnJ^ H + Ktf 

Xi Xj V a Xi l J l J 



where 



J 



jhH 

ij 

hH 

ij 



K, 



hH 

ij 



(s - 4m 2 ) 



%-y 



( 2 2 2 2 2 \ 

(s, m x ,m h ,m H , m^.m^o), 



,(4.9) 



(4.10) 



-4T l + 2(-2m 2 x + m 2 h + m 2 H )T -2y l 

n( A 2 \ ~\ 7 ~1 / 2 2 2 2 2 \ 

-2(s - 4m x )J^ j (s,m x ,m h ,m H ,m x o,m x o), 



(4.11) 



{-T 2 -(s + 2m 2 x -m 2 - m 2 H )T x - im\ - m 2 h )(m 2 x - m 2 H )T 
-y 2 + [-{ml + ml){m 2 x + m 2 H ) +4 m ^}( S , m >UL4,^)i 



(4.12) 
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Higgs (h, i7)-neutralino (x®) interference term: 



1=1 



E 

r=h,H 



(jhHr qXXT 

s — ml + iT r m r 



n xU h n x\}x H 



j — 2m x + [m x (m 2 H + m 2 h -2m 2 x - 2m^ ) 
+m x o (s - 4m x )]J 7 (s,m x ,m 2 H ,ml,m x o)^ . 



(4.13) 



As mentioned above, the couplings C s l 3 and C hHr (r = h,H), as well as all the other 
couplings appearing in this Section, are defined in Appendix A. The functions J-, 7*. and 
yk (k = 0, ... ,4), and all other auxiliary functions, are listed in Appendix B. By I\ and 
Th we denote the widths of h and H, respectively. 

The expressions for hh final state are obtained from the above by replacing mg, c hHr , 
Cg iXH , Cg jX with m h , C hhr , C^ iXh , C^ jX , respectively, and multiplying uj by a factor 
of 1/2 for identical particles in the final state. The contributions for HH final state are 
obtained in an analogous way. 



2. XX AA 



Similarly to the final state hH, this process proceeds via the s-channel CP-even Higgs 
boson (h and H) exchange and the t- and u-channel neutralino (x®, i = 1, . . . , 4) exchange 



~{h,H) . ~(x°) , ~(h,H- x °) 
WAA = W A A + A A + W AA : 



(4.14) 



• CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 
AA 



E °s 

s — 



r=h,H 



m, 



1 + i T r m r 



(s-Am 2 ); 



(4.15) 



neutralino (x®) exchange: 



w 



(x°) _ 1 
AA o 



where 



I- 



AA 



K 



AA 

ij 



T A A -j A A tt~ A A 

m v o m v oi 7 -« + m y m v oJ r ,' + A.--- 



, (4.16) 



(s-4raj) 



To-Eo 



/ 2 2 2 2 2 \ 

{s,m x ,m A ,m A ,rn x o,m x o), 



(4.17) 



= [4Ti + 4(mJ-m^)75 + 2yi + 2(s-4mJ)yo](s,m5,m^,mi,mJo,m^), 



{ - T 2 - [s + 2 (m£ - mi)]Ti - {m\ - m 2 A ) 2 % 

-y 2 + [4 m x - {m 2 x + m 2 A ) 2 ]y } (s, m% m 2 A , m^mjo, 



m x°< 



(4.18) 
(4.19) 
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Higgs (h, iZ)-neutralino (x®) interference term: 



w 



(h,H- x °) 



AA 



5> 



i=i 



(jAAr (jXX r 

. . . s — ml + iT r m r 

r=h,H 



(cf xA f 



: I — 2 m x + [2 m x {m\ — m 2 — m 2 p ) 
-m x o (s - Am x )]T(s, m x , m\, m\, m^o)}- 



(4.20) 



3. XX -*• fc-A 



This process proceeds via the s-channel Z and CP-odd Higgs boson (A) exchange as well 
as the t- and it-channel neutralino (x®, i = 1, ... ,4) exchange 



WhA - w hA + w hA + u> hj4 + u> hj4 + w hA + 
CP-odd Higgs-boson (A) exchange: 



(4.21) 



w 



(A) 
hA 



(jhAA (jXX A 



s - ml + i Fa ttla 



• Z-boson exchange: 



w 



(Z) 
hA 



(jhAZ (JXXZ 



s — m z + iTz mz 



m 7 s 



x jm^ 58 + s2 [6 fn x ( m A ~ m h) 2 ~ ^ m x m 'z ~ 2 m z ( m A + m h)] 
+sm 2 z [—\2m x (m\ — m\) 2 + 8m| to 2 {m 2 A + m^) 

+m| (m^ - m 2 ) 2 ] + 2m|m 2 (m A - m 2 ) 2 }; 



neutralino (x^ 1 ) exchange: 



(4.22) 



(4.23) 



(x°) 



where 



Z; 



J 



ij 
hA 



K hA 



4 

Y. c 



x°xh c x^xA f C $ xh y (Cp jXA 



s % - 3^o 
2(rr& 



/ 2222 2 \ 
{s,m x ,m h ,m A ,m x o,m x o) 



m h) 



%-y 



i 2222 2 \ 

(s,m ,m h ,m A ,m ,m 0j 



{ - T 2 + (2m 2 x + m 2 A + m 2 h - s)7i + 3^2 
+{m 2 (^ + m 2 - m 2 ) - m 2 A m 2 } [% - y ] } (s, m\ , m 2 , m\, m 2 ^ , m 2 o); 



(4.24) 

(4.25) 
(4.26) 
(4.27) 



- 10 - 



Higgs (A)—Z interference term: 



w 



{A-Z) __ 
hA 



2 Re 



QhAA (JXX A 



(jhAZ (jXXZ 



s — m\ + i Ta m,A j s — m? 7 + iTz mz 



m x (s- m|) (m\ - m 2 ) 



VTLr 



• Higgs (^l)-neutralino (Xi) interference term: 



(4.28) 



w 



(A-x°) 
hA 



i=i 



(jhAA (jXX A 

s - m 2 A + i Ta tua 



x |s m x q + m x {vr?A — mfy J r (s,m x ,m 2 l ,m A ,m' x o); (4.29) 



Z-neutralino (xf) interference term: 



w 



(z-x°) 
hA 



i=i 



s — m 7 + i Y z mz 



r,x°iXh n x°xA 



ITlr 



(s + 2 m\ -2m: — m\ - m 2 h ) m 2 z 



+2 



|m x o s [m 2 z m x o + m x {m\ — m 2 h )] — m 2 z m x m x o {m\ — m\) 

, 2 2 2 2/2, 2\ / 2 , 2\ , 2/2 2\2 

+m z m A m h - m z {m A + m h ) (m^ + m x ) + m z (m^o - m x ) 
+m 2 i (m 2 A -ml) 2 }j r (s,m 2 : ,ml,m 2 A ,m 2 x o) . (4.30) 

The expressions for HA final state are obtained by replacing mh, Cg iXh and C 1 ' 1 ^ in the 



y-yH 

above with mjf, C 5 4 and C MAZl , respectively. 



HAZ 



4. XX -> 

This process proceeds via the s-channel Z and CP-even Higgs boson (h, H) exchange as 
well as the t- and n-channel chargino , k = 1,2) exchange 



~(h,H) . ~{Z) 
WH+H- = W H+H~ +W H+H~ 



+ W 



(x±) 

H+H' 



+ W 



H+H- 



+ W 



(^-x±) . 

H+H- ■ 



(4.31) 



CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 
H+H- 



s — m 2 + iT r m r 

r=h,H 



(s-4m 2 x ); 



(4.32) 
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Z-boson exchange: 



w 



(Z) 

H+H- 



C 



H+H-Z (jXXZ 



s — m 2 7 + iTz mz 



chargino {x k ) exchange: 



(4.33) 



w 



(x±) 

H+H~ 



2 

E 

k,l=l 



where 



t h+h- 
l ki 



{{C H k C» + D H k D H t ) (-2m 2 x ) [% - y 

+ (C H k C H X - D H k D^)(s - 2m\) [% - y ] } (s, m 2 x ,m 2 H± , m 2 H± , 



(4.34) 



2 2 \ 

,m ±,m ±), 



3, 



H+H- 



kl 



rtH rtH 



- 4T ± - 4(m 2 x - m 2 H± )% 

23^i - 2(s - 4777^)3^0 {s,m 2 x ,m 2 H± ,m 2 H± ,m 2 x ±,m 2 x ±), 



K, 



H+H- 



kl 



(C? k + D* k D*i) -T 2 -(s-2 m' H± ) T x - (m x - m' H± f T 



2 \2 , 



(4.35) 
(4.36) 



+ (c?*c75-d*d3; 



+ 2m 2 (m x - m 2 H ±)y 



2m{T 1 -y 2 



■\-{m x -m 4 H± )y \(s,m x ,m% ± ,m 2 H± ,m 2 x ±,m 2 x ±), 

(4.37) 



and 



= \c** xH | 2 ± |Cp^ xH | 2 

D H k ^ c ^ xH ~{cf xH ~y± cf xH ~{cf xH ~y- 

Higgs (h, -£f)-chargino (xt) interference term: 



(4.38) 
(4.39) 



w 



(h,H- X ±) _ 
H+H- 



2^ Re 



k=i 



E 



(jH+H-r (jXXr 

s 

. s — m 2 + i T r m r 

r=h,H \ 



x <j Cf k (-2 m x ) + [C? k 2 m x (m 2 H± - m 2 - m 2 ±) 

+C H k m ± (s-4m 2 )] F(s,m 2 x ,m 2 H± ,m 2 H± ,m 2 ± ) \; (4.40) 
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• Z-chargino (xt) interference term: 



w 



H+H- 



k=i 



(jH+H-Z (jXX Z % 

s — m 2 z + iTz mz 



u +k 



x <j s - 2 m 2 H± -2m 2 x + 2 m 2 ± 



+2 



2i i 4 <-»/__ 2 i 2 \ 2 

sm i+m Hi - 2 (m v + m ±)mrr ± 
X k n *- x k n 

, / 2 2 \2l r/ 2 2 2 2 • 

+(m x -m x ±) J F(s,m x ,m H ±,m H± ,m x ± 



(4.41) 



5- XX -> W+ff- 

This process involves the s-channel CP-even (h, H) and odd (A) Higgs boson exchange as 
well as the t- and u-channel chargino (xt, k = 1, 2) exchange 



'(x±) 



CP-even Higgs-boson (/i, H) exchange: 



+ u; 



(A-X*) . 
WH 



(4.42) 



(h,H) 
WH 



E 

2 



to 2 + z r r m r 



4 m 



s 2 - 2 (m^± + m. 2 ^) s + (m^ ± - m 



2 \2 



m 



w 



CP-odd Higgs-boson (A) exchange: 



(4.43) 



w 



(A) 
WH 



C 



W-H+A (JXXA 



s — m? A + i Ta m,A 
s 2 -2 (m 2 H± + m 2 w ) s + (m 2 H± - 



nr 



2 \2 
W> . 



rn 



w 



chargino (xt) exchange: 



(4.44) 



w 



(x±) _ 

WH 



m 2 .4^. 



TO 4- TO 4- T" 1 *^ -4- TO TO 4- 7^ 4- K WH 
Xt xf H +m X m Y t J kl + 



x k 



(4.45) 



where 



l ki 



, n HW* n HW , n 
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jWH _ 
J kl — 



A M — 



% 



2 - (s - m 2 w - m 2 H± - 2 m 2 ) %_ + T + 6 m w m 2 x y 



+ (C^Cgr - D-rKi W ) 6 m^mi % - y 2 



HW* r^HW\ 



(s,m ,m H± ,m w ,m ±,m ±), 



(4.46) 



Re(C^ k w *C^ w + D 



2T 2 + 2(2m2 -m^)Ti 



-2(m 2 - rr^)(rr? x + 2 m^)7^ - 3 m 2 ^ ^i - 
+ Re{C^C^ - D™*D™) [6m^r 1 + 6 m 2 w (m 2 - m 2 ^ )T 
-2J 2 -(s-mH±-2m^i 



r.y(3) v 



(22 2 2 2 \ 

(s, m x , m H± , m w ,m v ± , m v ± J, 



X; 



(C W* C W + D HW* D HW } \ {s _ m 2 x _ 2m 2 w)T2 _ G T(2) ^ 



- G^g T + ( s - 2 m^) 3^2 + G' W ( J> y 
+ (G^CT " D™*D™) \ 6 m^m 2 7i - m 2 y 2 



y(4) 



(4.47) 



and 



G WH = ( m x ~~ m H±) ( m x ~ m w) ( m X + 2m w), 



Lx 
Lx 



Y(l) 

y( 2 ) 

Y(3) 
WH 
Y(4) 



G 
and 



Y(5) 
WH 



/ 2 2 2 2 2 > 

(s, m ,m H ± , m w , m ± , m ± , 



Xfc X; 



G^ = s(m 2 + 2771^) -m x - m 2 x (m 2 H± + Sm^r) - m 2 v m 2 H± , 
G WH = s ( m x + 2r77 w) ~ 2m i + m 2 (m 2 H ± - mw) - 2m 2 v (m 2 H ± + m^), 



x ' x v 

2 2 \ i 2 



G 



T(l) 



3 s m 2 ^ - 12 m^m 2 + 3 m^ra 2 ^ - 3 m^, 



(4.48) 



s 2 - s(2m 2 + 2m^± + 3m^,) + 2m x + (2m 2 + m 2 H± )(m 2 H± + 3m^) - 2m^, 

s (m 2 - 2m 2 v )(m 2 x - m 2 H ±) + 4m 2 v m x 

+ m 2 x (m A H ± - 5m^m^± + 2m^) - 2m w m 2 H ±, 

sm 2 (m 2 - m 2 ^) - m x — m x (m 2 H ± + 3m^) + m 2 x m 2 v (2m 2 H ± + 3m 2 ,/), 



x^" 



C 



± c. 



xtxH- ( c xtxW 



V 



(4.49) 
(4.50) 



Higgs (h, i?)-chargino (xt) interference term: 



w 



WH 



-2^E E 



fc=l r =h,H 



s_ 

s — ml + i T r m r 



m 



W 
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riHW* tt(+) wh _l r<HW* tt(-)wh 



(4.51) 



where 



2 m x m ± (s + m 2 ^ — m^± ) + m x m ± s 2 + 2 s (m ± — m 2 — m^i 



,4 4,o/'2 2w2 2 \ -r-/ z ^ s , 

+m H± - + 2 (m ± - m x ) (m^ -m H± ) F(s,m x ,m H± ,m w ,m v ±), 



x 

2 _2 



2,/2, 2 . o„2\ , o 2 /■ 2 2 \ 

-s + s(m w +m H ± +2m x ) +2m x (m w -m H ±) 

+ | - s 2 m 2 ± + s[2 m 2 (™< 2 ± - m 2 ) + m 2 (m 2 H± + 3 m 2 ^ ) 



+m ±(m H± +m w ) - 2m w m H± ] 

Xk 



+m 2 (mly - m 2 H± ) ( 2m - x ± ~ 2m x + m ^ ± ~~ 3m w)} 



x^(s, m 2 m^i , m 2 ^, m 2 ±); 



2 

(4.52) 



(4.53) 



Higgs (^4)-chargino (x k ) interference term: 



where 



WH 



-2ReH 

k=l \ 



(jW'H+A (JXX A 
s — m? A + i Ta rriA 



m 



w 



[D^A { k +)WH + D™* A (-)WH 



(4.54) 



a(+)WH 



A 



(-)WH 



2 s (m 2 H± + m 



2 \ 
w) 



,/2 2 \2 2 2 2 2n 

+(777,^-777,^) F(s,m x ,m H ±,m w ,m±), 
s (s - m 2 ^ 



(4.55) 



2m 2 v m 2 H± ] +m 2 c (m 2 H± - m 2 v ) 2 ^J r {s, m 2 , m 2 H± , m^, m x ±)- (4.56) 



Note that wwh in eq.( 4.42j ) does not include the contribution for W~H + final state. The 
contributions to the W~H + and W + H~ final states are obviously identical and so the 
total contribution for W^H^ is twice that of W + H~ . 



6. XX Z/i 



This process proceeds via the s-channel Z and CP-odd Higgs boson (A) exchange as well 
as the t- and u-channel neutralino (x®, i = 1, . . . , 4) exchange 



WZh 



~(A) . ~(Z) . ~(x°) , ~{A-Z) 
W Zh + W Zh + W Zh + W Zh + w ™ + W 



Zh w Zh 



(4.57) 
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• CP-odd Higgs-boson (A) exchange: 



w 



(A) __ 
zh — 



(jZhA (jXX A 



s — m A + i Ta rriA 



s 2 — 2 (ml + m|) s + (m| 



m 



2 \2 

z) . 



rrir 



Z-boson exchange: 



w 



(Z) 
Zh 



12 



(jZZh (jXXZ 



s — m 2 7 + i Tz mz 



m z s 

,2 \ „3 



{m 2 x [6 s 4 - 12 (m\ + 2 m|) s 3 + (32 m% + 12 m 2 z m 2 h + Qm\) s 2 
- (64 m% - 20 m| m£ + 12 mfj sm 2 z + 2(m 2 h - m 2 z ) 2 m\\ 
+m|[s 3 + (10 m| - 2m 2 h ) s 2 + (m| - m|) 2 s]} ; (4.59) 



neutralino (%°) exchange: 



(x°) 



mi ' 



xh h n x °i xh * 
s b s 



m v p m v ql? h + m x m x oJ^ h + 



where 



(4.58) 



(4.60) 



- T 2 - (a - m| - m| - 2mJ)Ti + G 1 ™ % + 6m 2 z m 2 x y 
HC^C** - ClCl*)\6m 2 z m 2 x T - y 2 



I 2 2 2 2 2 \ 

(s, m , m h , m z ,m o,m o), 



J, 



Zh 



(C$ t C** + C\f%) \-2T 2 + 2(2m 2 x - m|) 7] 



+(C&C$*. - C*Cg) [ 6 m| Tx + 6 m 2 z (m 2 x - m 2 )% 

- 2y 2 -(s-m 2 -2m 2 z )y 1 



-G Y Zh 3) yo 



I 2 2 2 2 2 \ 

(s,m x ,m h ,m z ,m x o,m x0j 



Kfj h = j (C* C** + C* C%) [ ( a - m\ - 2 m|) T 2 - G T z f T x 



- G T z f% + (s-2m 2 z )y 2 + G Y z ^y, 



HC^) ~ C\fi%) \ 6 m|m 2 T x - m\ y 2 



+g y z ( h y° 



I 2 2 2 2 2 \ 

(s,m x ,m h ,m z ,m x o,m Yih 



(4.61) 



2(m 2 x - m 2 z )(m x + 2m z )% -?>m z y x - G Y z f y 



(4.62) 



(4.63) 
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and 



and 



G 



C 1 y z > — 

Livi — 



G 
G 
G 
G 
G 



T(l) 
Zh 
T(2) 
Zh 
T(3) 
Zh 

y(i) 

Zh 

y( 2 ) 

Zh 
Y(3) 
Zh 
Y(4) 
Zh 



s(m 2 + 2m|) - m x — m 2 x (m 2 h + 3m|) - m z m 2 h , 
s(m 2 x + 2m|) - 2m x + m 2 x (m 2 h - m|) - 2m 2 z {m\ + m|), 

/ 2 2\ / 2 

(m x - m J (m x 



"x 1 x v 
m|) (m 2 + 2m|), 



G 



T(l) 
Zh ' 



3 s mfz — 12 m 2 z m 2 + 3 m^m| 



3m ; 



s 2 - s(2m 2 + 2m 2 , + 3m|) + 2m x + (2m 2 + m 2 h )(m 2 h + 3m|) - 2m 
s (m 2 — 2m|)(m^ — m|) + 4m 2 z m x 



+ "M m h 



hm 2 z m\ + 2m^) 



2m z m 2 h , 



G Y Z ^ = s m 2 (m 2 — m 2 z ) — m® — m x {rr? h + 3m|) + m 2 m z (2m 2 l + 3m|), 



r x - ^iX-x^" 



x°x^ 



• Higgs (A)-Z interference term: 



71' 



(A-Z) 
Zh 



Re 



(jZhA (jXX A 



(jZZh 



s — m 2 A + i Ta rriA J s — m 2 z + iT z m z 
m x (m| - s) 2 



[s 2 -2(mi + m 2 z ) S + (m 2 h -m 2 z Y]; 



Higgs (A)-neutralino (Xi) interference term: 



(A-x°) 
Zh 



2 /Z Re 



i=i 



(jZhA (jXX A 



b S °A 



ITLr 



s - m 2 A + i Ta vtia 
S 2 + (m 2 + m 2 ) S -{ S 2 m x o(m x + m x o) 
+S [2m 2 m|-(m 2 +m|)(m x + m x o) 2 ] 
+m x (m x + m x o) (m 2 - m|) 2 } F(s, m 2 , m 2 , m| , m 2 ? ) 

Z-neutralino (x° ) interference term: 



(z-x°) 

Zh 



i=i 



^ZZh £<XXZ 

s — m| + iT z m z 



) 



siXiXh n X°xZ 



m? 



x jm x [-2s 2 + s(2mj + 3m|) + m| (2m 2 - 9m| - m 2 - 2m 2 )] 



+77V m\\-s + {2m 2 +m 2 z + m 2 -2m\)\ 



m x {s 2 m 2 x o + s [-{m 2 x + m£o) (m| + m|) + 2 m| m|] 

-2m| + 5 m 2 , m z + 4m^ m| + m| (m?. - m^ ) 2 
-m| m| (3 m| + 2 m 2 ,) + m 2 , m^} 

r O O / 00 00 /1\ 

+m y p [s m y + s (— 3m y m z + m z m x o — 2 m x m ft — 2m^J 



x 

+10 m 2 x m z - m x o m% + m z (m 



x ' 

,2 4 



2 \2 



22/2 2 2 \ , 2 Ai -t-( 2 2 2 2 

+m h m z (m z - m - m ) + m m ft J ^(s, m , m h , m z , m 



(4.68) 



The expressions for Zff final state are found by replacing h with H in the above. 



7. XX — ZA 



This process involves the s-channel CP-even Higgs boson (h and H) exchange and the t- 
and ^-channel neutralino (%°, i = 1, . . . , 4) exchange 



w Z A = w ZA ' + w ZA > + w ZA 



-Xx°) 



(h,H- x °) 



(4.69) 



• CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 
ZA 



(jZAr C^ XV 



s — m 2 . + i T r m r 
s 2 - 2 (m 2 A + m|) s + (m^ — m|) 2 



r=h,H 
2 



(s-4m 



771: 



(4.70) 



neutralino (xf) exchange: 



w 



(x°) 

ZA 



1 V- C X?XA C X°X^4* 



7/4 7/4 7 A 

m v o ra^o 7 ? - „• + m v m v o J 7 a + „■ 

Aj Xj l J A Xi 4 J L J 



(4.71) 



where 



- T 2 -{s-m 2 z - m\ -2m 2 x ) T x + T + 6m 2 z m 2 x y 
HC^C** - C\£%) [ - 6 m 2 z m 2 x T + y 2 

r^Yil) 2 2 2 2 2\ 

~ G ZA SO [{ s , m x> m Ai m Z, m x o, m x o), 



(4.72) 
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and 



J 



ZA 



(C* C** + C\fl%) [2T 2 -2 (2 - mi) 7i 

+2(m 2 - m|)(mj + 2m|)T + 3m 2 z y 1 + G^ 2) 



+(G*,G** - C**) 6m|T 1 + 6m|. 



J Ai^ Aj 



m. 



m A )T Q 



-2^ 2 -(s-m|-2m|)^i 



^(3) v 
G ZA ^0 



(s,mj,m^m|,ra\,m^), (4.73) 



if, 



ZA 



l U Vi U Vj Ai A 



Ai^AjJ 



G T z f T + (a - 2 m|) y 2 + G^ % 



2m|)T 2 -G^ 2) T 1 



y(4) 



+(C^G** - G*-G**) -6m|m^r 1 + m^2 



GS ) = «K + 2m|) 



m. 



G 
G 
G 
G 
G 
G 



T(2) 
ZA 
T(3) 
ZA 

y(i) 

ZA 
Y(2) 
ZA 
Y(3) 
ZA 
Y(A) 
ZA 



X 

s(m x + 2m|) — 2m x + m 2 (m/4 — m z ) — 2m z (m A + m z ), 
(m 2 - m\) (m 2 - m|) (m 2 + 2 m|), 



m 2 (m 2 A + 3m|) 



(s, m x ,m A , m z , m^, m^o), 



(4.74) 



2 2 

m z m A , 



G 



T(l) 
ZA ' 



G 



Y(5) 
ZA 



3 s m?7 — 12 m z m x + 3 m z m A — 3 m-i, 

s 2 - s(2m 2 + 2m A + 3m|) + 2m^ + (2m 2 + m 2 A )(m 2 A + 3m|) - 2m|, 

s (m 2 — 2m|)(m 2 — m A ) + Am 2 z m x 

+ m x (m A — hm 2 z m 2 A + 2mi) — 2m z m A , 

s m 2 (m 2 — m 2 z ) — m x — m x (m A + 3m|) + m 2 m 2 z (2m 2 A + 3m 2 z ). 



"X^ X A> X ""XV A ' "■■"Z,/ ' ■•"x 

G T ZA ^ and G Y ZA ^ are obtained from G^fc 3 ^ anc ^ ^zfc 1 5 ^ by replacing m^ with m^. 
Gyj and C Ai have already been defined for the Zh final state in eqs. (4.64) and ( |4.65| ). 

• Higgs (h, i^)~neutralino (xf) interference term: 



w 



(h,H- x °) 
ZA 



-2^ Re 



i=l 



E 



(jZAr fjXXr 

S 



r=h,H 
2 i ™2 



r,x°x z ^.x°xA 
°A °P 



1 



s — m 2 r + i T r m r 
s 2 + (m A + mi) s + 2 m x (m x o — m x ) (— s — mi + m A ) 

+ | — (m^o + m x ) (s — m^)[m x o s + 2m x (m x o — m x ) 2 — m x m^] 
+mi m x [—2 (m x o — m x ) 2 (m x o + m x ) + 4 m x m\ 

+m 2 z (m o — 3 m x ) + 3 s m x ] + s mi (m 2 — 2 m 2 ^) > 

'M Xi ) 
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xJ 7 (s, m 2 m\, m 2 z , m 2 ) 



(4.75) 



8- XX -> WW 



This process involves the s-channel CP-even Higgs boson (h and H) and Z exchange, and 
the t- and n-channel chargino (x^, A; = 1,2) exchange 



~{h,H) . ~{Z) . ~( x ±) . -(KH-x^) . -(Z-x*) 

www = w ww ' + w\yfo + wfttf + w ww + w ww A : 



(4.76) 



CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) _ 
WW — 



E 

r=h,H 



(jWWr C^ Xr 



s — ml + i T r m r 



(s-4«v)- 



2 - 4m w s + l2m 4: 



i m 



w 



(4.77) 



Z-boson exchange: 



w 



(Z) 
WW 



(jWWZ (JXX Z 



s — m z + iTz mz 
s 3 + lQm w s 2 



(s — Am 2 



68 m, 



48 



12 



(4.78) 



chargino ) exchange: 



w 



(x±) _ 
WW 



4-E 



fc,/=i 



rWW , , t/WW 



(4.79) 



where 



A kl 



T ww 
J kl 



(s-4 rr&) T 2 - G^ w Ti + T + (a - 4 m 2 w ) ^ - 3> 



+(C^*C^ - D £D 



I8m 4 w ml% 



-ml(-s 2 + Am 2 w s + Q m* w ) y \ (s, m 2 x ,m 2 w , m 2 w , m\ ,m 2 ± ), 

Xk Xi 



/ ^iW*/~<W 



12 T 2 + 12 m 2 ^ (m 2 ^ - 2 m 2 ) Ti 



(4.80) 



+I2m w (m x + m^rw 2 . - 2m w )% - 4 (s - m w )y 2 

+2m 2 v (2s -3m 2 v )y 1 + G^^y Ms,m^,m^,m^,m^±,m^±), (4.81) 



- 20 - 



K, 



WW 

ki 



(C%C% + D™D%) [-%-(s-2m w -4 m 2 ) % - T 2 
~G\vw r i " ( m w ~ ™') 2 ( 2 m w + ™\f % + m 2 x (s-4 m 2 w ) y 2 - gQ w y 



+{C^CZ ~ D%D%) - 18 m 4 w m 2 x %_+ y 4 + (2m x -4s m 2 w + m 4 w ) y 2 



-G 



Y(4) 
WW 



3^o 



[s,m x , m w , m w , m ± , m ± ) , 



(4.82) 



and 



and 



G ww = 2 i s ( m x + 2m w) -44K + 4)1. 

= s(m 2 + 2m w ) 2 + m w (-4m x - 10 m 2 -Amfy), 



G 
G 
G 



T(3) _ 

~~ 

T(4) 
WW 
Y(l) 
WW 



—2s (m 2 x + 2 m 2 ^) + 6m*+6m^ m 2 ^ + 5 ra w , 



! ra w -6mJ m 2 ^ 



4m°, 



sm 2 (m 2 + 2™^) + 6 m 2 mfy + 4m w (m 2 - m w ) 2 , 



G## = 2[2rr& s 2 + s(2m 2 - 2 - 11 m^) 
+2 m 2 w {m x + 4m 2 mfy + m^-)], 

G ww = m x( m w ~ m x )[ s ( m x + 3m w) - 2ra w( 2m x + m w)]' 

G ww = (m 2 w - m 2 x )[s (8^ -4m 2 m 2 ,) - m 2 (8^ - m 4 x - m 2 m 2 w) ], 



G 



w 

±k 



Higgs (/i, i?)-chargino (xt) interference term: 



± cf xW 



(eg**-)'; 



(4.83) 
(4.84) 



where 



H, 



(+)ww 



w 



WW 



k=l 



< + c w k Hir )ww 



qWWt C^ XT 

Tr s — ml + iT r m r 

,r=h,H 



E 



m 



w 



H (-)ww = 



-m ± (s 2 - 2 m 2 ^ s) + m ± 
+3 m 2 ^ m 2 ) - 24 m 2 



(4.85) 



^k ^k 

+2 m y \ s 2 m 2 ± + s \2 mtxr — 3 m\j m 2 — m 2 ± m 2 ^ + (m 2 — m 2 ± ) 2 1 

I x fe * x fc * x fe 

-4 m w + 2 m w (m 2 + m 2 ± ) 

Xfc 

+2 771^ (m 2 - m 2 ± ) 2 \T{s,m 2 x ,m 2 v ,m 2 v ,m 2 ± ), 

- (m 2 + m 2 ±) s 2 + 2 s (2 m™, + mlm 2 ± 
A x fc x fe 

T(s,vr?, m 2 w ,m 2 w ,m 2 ± )\ (4.87) 
A x fc 



(4.86) 
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Z-chargino (xt) interference term: 



w 



WW 



1 



2 



' qWWZ (JXX z 

s — m? z + iTz mz 



n w 
u +k 



x< s 3 + s 2 (18m w - m 2 - 3m 2 ±) 

f2s [— 14m^/ — 18m 2 ni^y + 6ra 2 ± rri^y — 3 (m 

2 \2 



Xk 



— m w (11m 2 — 3m 2 ±) — 6m^] 

X k 

s 2 m 2 ± (m 2 — m 2 ± + 4m.w) + s[8m w — 5m w (m 2 ± + 3m 2 ) 

1 Xk " X k A 



+Am 2 v [3 (m 2 - m 2 ± 

Xk 



+6 



L x fe v x Xk vv ' L 

+2m 2 w (3mi-5mlm 2 ± +2m 4 ± ) + (ml - m 2 ± ) 

x x x k x k x x k 

+m w [Am w -2m w (3m 2 ± + 5m 2 ) 

Xk 

i o 2 2/2 2 \ i o I 2 2\3ll <rv 2 2 2 2 \ 

+8m w m x (m x - m^± ) + 2 (m^± - m x J JJ.F(s,m x , m w ,m w ,m^±) 



(4.88) 



9. XX — ZZ 



This process involves the s-channel CP-even Higgs boson (h and H) exchange and the t- 
and ^-channel neutralino (%°, i = 1, . . . , 4) exchange 



. ~(x°) , ~(h,H- x °) 

w Z z = w Z z + w zz + w zz ■ 



• CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 

zz 



E 

r=h,H 



(jZZr QXXr 

S 



s — m 2 , + i T r m r 



(s — Am 



2 s 2 - 4m| s + 12 m\ 



16 m\ 



neutralino (x®) exchange: 



(4.89) 



(4.90) 



where 



tZZ 

ij 



w 



(x°) 

ZZ 



2 m 



»J=1 



m v o m^l zz + m y m^J zz + .FQ? Z 

Xj t J ^ Xi l J L J 



(4.91) 



x '(a - 4 m|) T 2 - Ti + G^ 2) T + (a - 4 m 2 ) y 2 - G Y z f y 
+ {C z *C Z j - D z _*D z _j) \ - \8m\m\ % 



-22- 



and 



K 



zz 

ij 



-m 2 (-s 2 + 4m| s + 6m%)y >(s,m 2 m 2 z ,m 2 z ,m 2 o,m 2 o), 



CZ* fiZ 
-i 



12 m| T 2 + 12 m| (m| - 2 m 2 ) 71 



+12 m| (mj + m 2 m| - 2 m|)T - 4 (a - m|)^ 2 
+2m| (2s - 3m|)^i + G^f ) D^o (s,m^m|,m|,m^,m^o), 



^ZZ J 2 



T 4 -(s-2m|-4m 2 )T 3 
-(m| - ) 2 (2 m| + m 2 f % + y 4 



+G^l 3) 3^2 - G^p 1 3^0 



Y(4) 



(s, m x , m z , m z ,m x o,m x o), 



(4.92) 



(4.93) 



(4.94) 



G 
G 
G 
G 
G 
G 

G 
G 



T(l) 

T(2) 
ZZ 

T(3) 
ZZ 
T(4) 
ZZ 

y(i) 

ZZ 

r( 2 ) 

ZZ 



y(3) 

ZZ 
Y(4) 
ZZ 



2[s(m 2 + 2m|) - 4m|(m 2 + m|)], 



= s(m 2 + 2 m|) 2 + m 2 z {-\m\ 



10 m 2 m 2 z 



4m 4 z ), 



-2 s (m 2 + 2m|) + 6mJ + Qm 2 m 2 z + 5m z , 



2 \2 



s (m + 2m z ) 



8m| + 18 m% m 2 



c 4 2 



-s m. 



m +2m z )4- 6m m z + 4ra z (m 



4 



4 m 



2 \2 

™z) . 



2 [2 m| s 2 + s (2 m 2 x m 2 z -2m x - 11 m| 

+2 m| (m x + 4 m 2 m 2 z + mf )] , 
s(m 2 x - 4m|) + 2m x - 4m 2 m| + m|, 
(m|-m 2 )[s(8m|+m 4 ™ 2 



m v m z) + Tn 2 (—10m z + m x — 3m 2 m 2 z )], 



and 



Cf f EE |C* Z | 2 ± |C-f lZ | 2 , 



D f jS C* z (c* z )- ± c5^(c* z )-; 
Higgs (/i, .ff )-neutralino interference term: 



(4.95) 
(4.96) 



where 



w 



{h,H- x {) ) 
ZZ 



5> 

i=i 



s — 



s — m 2 . + i T r m r 

r=h,H 



2m| 



(4.97) 



= m x {s 2 -2s(m 2 + m| - m 2 o) - 4m| (m 2 -2m|-m 2 o) 
+2 [s 2 m 2 o + s {(m 2 - m 2 o) 2 - m| (3 m 2 - 2 m| + m 2 o)} 
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+2m| {(m 2 - m 2 ) 2 + m| (m 2 - 2m| + m 2 )} 

Aj A Aj 



(4.98) 



H 



(-)ZZ 



-m 



{2 2 
s — 2m z 



xJ(s,m^m|,m|,m^)}, 

s 2 (m 2 + m 2 ) - 2 m| (2 m| + 3 m 2 + m 2 ) s 

+24m|m 2 F(s, m\ , m|, m|, m£o)}. (4-99) 



Note that = 0, because Cy iXj is pure imaginary and C^ 1 * 3 is real. 



10. XX — // 



This process involves the s-channel Higgs boson (h, H and A) and Z boson exchange and 
the t- and u-channel sfermion (f a ) exchange 



(h,H) . ~(A) . ~(Z) . -(/) . ~(h,H-f) . ~(A-Z) . ~U-f) . ~(Z-f) 

If + w If + w If + w If + 4/ + 4/ + w If + 4/ : 



CP-even Higgs-boson (h, H) exchange: 



w 



(h,H) 

Is 



r=h,H S 



s — m 2 . + i T r m r 



(s -4m^) (s -4m 2 ); 



• CP-odd Higgs-boson (A) exchange: 



w 



(A) 
If 



r<ff A nxx& 

Op Op 



s - m\ + i mA 



• Z-boson exchange: 



w 



(Z) 
If 



s — m 2 z + iTz mz 



Jl2\C{f Z \ 2 ^(s-m 2 z ) 
I m z 

\C{/ Z \ 2 (s + 2m}) + \C f A fZ \ 2 (s -Am 2 f )\ (s - 4m 



,2 \2 



+ 



• sfermion (/„) exchange: 



(4.100) 



(4.101) 



(4.102) 



(4.103) 



w 



(/) 
If 



-Y, 

a,b 



[C a + C b + + D%D\) [r 2 - 2 (m 2 + m 2 ) 7} 

+ {(m 2 +m^) 2 + 4m 2 m}}r 
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+{CtC b _ - D a _D b _) [T 2 -2(m 2 x + mj) T x + (m* - m 2 ) 2 T 
-{ClC h + - D a _Dt)4m 2 x m 2 f y 

+{c a +c h + + cic b _)sm 2 x y + {c a + c b + - c a _c h _) s m) y 

+ (D a + D b + - D a _D b _) [y 2 + (m 2 , + m 2 f ) 2 y 
+{C a +D b + + D a + C b + )m x m f U {m 2 + m 2 )T - 471 + s y 



(C%D b + - D a + C b + )m x m f y l 



I 2 2 2 2 2 \ 

(s,m x ,m f ,m f ,mj^mjj, 



(4.104) 



where 



c a ± 



I a a 1 2 _i_ I a a 1 2 
\ A fL\ ± \ 1V fR\ i 



= A a /L (A^)* ± (A? L )*A^. 



(4.105) 
(4.106) 

6 for squarks 



In the above, a is the index for sfermion mass eigenstates so that a = 1, 
and charged sleptons, and a = 1, 2, 3 for sneutrinos. 

The symbol / represents each fermion: f = u, c, t, ■ ■ ■. The symbol / should be 
understood as follows. For up-type quarks, down-type quarks and charged leptons, the 
corresponding symbol f a represents u a , d a and e a (o = 1, • • •, 6), respectively. For neutrinos, 
f a represents v a (a = 1, 2, 3). For eg. for 4?, the last argument at the end of eq. (|4.104| ) 
should read (s, m x , m 2 ,m 2 ,m 2 la ,rn'~ b ) when a = 1, • • •, 6. 



Aa _ \(u)L Aa 



The coupling A°j L for each fermion is defined by A^ ai in Appendix A as A^ L = A^^, 

^3ai ' i e * c - Similarly, the 



<u)L 



A (n)L A a - 

i4 -3al ' ly dL 
(f)R 



Il lal > iV sL 



A (d)L A a 



coupling Ay R is defined by A, 



Higgs (h, H )-sfermion (/ a ) interference term: 



w 



(h,H-f) 
If 



E 



=h,H 



s — m 2 + iT r m r 



C\2m x rrif 



{2+ [s-2( 



m x + m 



f -mj)] 



+D a + {s+[s( 



^T(s,m 2 x ,m 2 ,m 2 pm~)^ 



2 2 2 > 



xj r (s,m x ,m 2 ,m 2 ,m 2 jj^ 



(4.107) 



Higgs (A)—Z interference term: 



'?/; 



(A-Z) 



Re 



r<ff A nxxA 

Up Op 



r<ff z nxxz 

U A °A 



s — m A + i r^4 rriA I \ s — + i 



X ^L S ( m |- S ); 



(4.108) 
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Higgs (j4)-sfermion interference term: 



w 



(A-f) 
If 



r ffA r xxA 

y^ Re C p C p 

^— ' s — m a + il a rriA 

a L Ji 

C+2 m x nifsJ : {s, m x , mj , mj , m' 
+£>+{ - s + s [m 2 x + mj -m~) 

xT{s,m x ,mj,mj,m 2 ~ ) j 



(4.109) 



Z-sfermion (/„) interference term: 



w 



(z-f) 
If 



r<ff z nxxz 
s — m% + iTz mz 



Zi + W + Re 



s — m\ + iTz mz 



z (-)ff 



(4.110) 



where 



= Cl{s + 2{m 2 x + mj-m 2 j)- - 2 s (mj, + mj) + 2 (mj + mj - m|- ) 2 



+8 m 2 mj 



-EY 2 2 2 2 \1 

+ (i " 3 ) (ml + m}-m 2 j)\F(s,ml,m},m 2 f ,mjj} : 



(4.111) 



2 2 

i f -mr 



+ {2s {m 2 x - mj) - 2 (m 2 x + mj - m~ ) 2 + 8m x mj] 
xj : (s,m x ,mj,mj,m 2 j )j. 



(4.112) 



This completes the list of all the tree-level two-body neutralino pair-annihilation chan- 
nels in the MSSM. 



5. Partial Wave Expansion 

In the literature one still often uses the usual approximation in terms of the expansion in 
powers of x (or, equivalently, WIMP velocity-square), (ctvmoi) — a + bx. This is because 
in the early days it was often much easier to compute the coefficients a and b, rather 
than the cross section itself [p^j . Furthermore, the partial wave expansion gives a rather 
good approximation to the exact result but only far enough from s-channel resonances and 
thresholds for new final states. (For a recent detailed study, see Ref. [p7gf| . ) 



- 26 - 



The expansion of (crt! M0 i), as defined in eq. (|2.2|) , is given by [gl| 16] 



(o"«M0i) = — T 



it) 



- (2w - w') x + 0(i 



a + bx + 0(x 2 ) 



(5.1) 



where u/(s) denotes dw(s)/d(s/4m 2 



Analogously to the function w(s) (eq. ( |4.3| )), the coefficients a and b need to be summed 

(5.2) 



over all possible final states /i/2- One can write them as [21] 



c6 {Am 2 , - (m h + m h f) v hh a flh , 



hh 



Y c6 (Am 2 , - (m fl + m h ) 2 ) v hh I b hh + 

hh y 



(5.3) 



-3 + -v f f 

4 71/2 



/ m /i+ m / 2 
2m| 



+ 



TO 



2 N2 % 



8m x 



where the summation extends over all possible two-body final states /1/2, the coefficient c 
is defined in eq. (4.4), and 



v hh = f3f(^m 2 x ,m fl ,m h ), 



(5.4) 



and the velocity Pf(s,mf 1 ,mf 2 ) was defined in eq. fl4.5|) . The "reduced" coefficients Oft/a 
and fe^yjj are given by 



a hh 
b 



327TTO?, 

A 



/l/2 64ttto2 



(5.5) 
(5.6) 



where u;(s) was defined in eq. (4.6) and w'(s) = dw(s)/d(s/Am 2 ). 

In this Section, we provide a set of expressions for the coefficients a and b in the case 
of equal- mass final states. Using eq. ( |5.lD , we have derived the coefficients for all the final 
states by using the analytic expressions for w(s) presented in the previous Section. In the 
case of unequal masses of the final state particles the resulting formulae are exceedingly 
lengthly and we will not include them here. 

In the literature one can find several analytic formulae for the expansion coefficients, 
including [12, [l3|, [U], |2|, but, due to different conventions and complexity of expressions, 
comparison is not always doable. We have checked our results for the a-coefficients in 
appropriate limits against published results and agreed in all cases. 

All the couplings and auxiliary functions appearing below are listed in the Appendices. 
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1- XX -» hh 

This process involves the s-channel CP-even Higgs boson (h and H) exchange and the t- 
and u-channel neutralino (x°,i = 1, . . . , 4) exchange 



CP-even Higgs-boson (h, H) exchange: 



a 



(h,H) _ 

hh — 

u(h,H) __ 

J hh — 



o, 

3 
64T 



E 



(jhhr (jXXr 



4 m£ — m 2 . + i T r m r 

r=h,H * 



neutralino (xf) exchange: 



r(x°) _ 



hh 



r(x°) 

J hh 



o, 



L E(^ X ") 2 (^*) 2 



16 7T 



x [4 m 2 (m 2 - ml) 2 + 4m x (m 2 - mfj (m x + m x o) A 
+3 (m x + m x o) (m x + m x o) A^ A/y 
where A^ = m? h — m 2 , — m^ . 
• Higgs (h, ii~)-neutralino (xf) interference term: 



fii 



(h,H- X °) _ n 
h hh 



hh 16 TT^ 



E 



(jhhr (jXX r 



. 4 ml — m 2 , +iT r m r 

r=h,H * 



c 



XiXh^XiXh 

s °s 



[2m x (ml - ml) + 3 (m x + m x o) A hi ] 



A li 



(5.7) 
(5.8) 



(5.9) 
(5.10) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



The expressions for HH final state are obtained from the above by replacing m^, C hhr , 
Cf x \ Cf xh with m H , C HHr , cf xH , C$ XH , respectively. 



2. XX -> AA 



Similarly to the final state hh, this process proceeds via the s-channel CP-even Higgs boson 
(h and H) exchange and the t- and w-channel neutralino (x* j* = 1> • • • >4) exchange 



7?,, _ ~(^) , ~(x°) , ~(^-x°) 



bAA 



L(h,H) 
AA 



(x°) , T(fc,H- x °) . 



+ &aY + & 



(5.15) 
(5.16) 
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CP-even Higgs-boson (h, H) exchange: 



a 



(h,H) 
AA 

Z(h,H) 
AA 



neutralino (xf) exchange: 



0, 

3 
64^ 



E 



qAAt qXX t 



4 m 2 — ml + i T r m r 

r=h,H X 



(5.17) 
(5.18) 



r(x°) _ 



l AA 



r(x°) 

AA 



o, 



1 v (r x ° xA \ 2 (r^ xA *\ 2 1 



16 7T 



«J=1 



A 2 A 2 
^Ai ^Aj 



x |^4 m 2 (m 2 — m/4) 2 + 4m x (m 2 — m 2 ^) (m x — m x o) Aa% 
+3 (m x - m x o) (m x - m x o) A Ai A A j , 



where A A i = ml A — m 2 — m 2 



Higgs (h, .ff )-neutralino (xf) interference term: 



(h,H- x °) _ n 
A AA u - 



r{h,H- x °) _ 1 



J AA 



16 7T 



i=i 



E 



(jAAr Cg Xr 



. 4 m 2 — m 2 + i T r m r 

r=h,H * 



[-2 m x (m? x - m 2 A ) - 3 (m x - m x q) Am 

A 2 
^Ai 



(5.19) 



(5.20) 



(5.21) 



(5.22) 



3. XX — 



This process proceeds via the s-channel Z and CP-even Higgs boson (h, H) exchange as 
well as the t- and n-channel chargino (x^, k = 1,2) exchange 



a H+H- 



~(h,H) ~(Z) ~( x ±) ~{h,H- X ±) ~{Z- X ±) 

a H+H- + a H+H- + a H+H- + a H+H- ! " 



b 



(h,H) 
H+H- 



+ ) + X h ' H - x ) +b {z ~ x±) ■ 



H+H- ' 
H+H- 



CP-even Higgs-boson (h, H) exchange: 



(5.23) 
(5.24) 



l H+H- 



U H+H- 



= o, 



3 
32T 



E 



C 



H+H- 



c xxr 



tt 4 to 2 — ml + i F r m r 

r=h,H 



(5.25) 
(5.26) 
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• Z-boson exchange: 



a {Z) - 
a H+H- ~ 

°H+H~ ~ 

chargino (xt) exchange: 



0: 

1 

4-7T 



C 



H+H~Z (jXXZ 



4 m| — m 2 z + iTz mz 



(m 2 - m 2 H± ); 



(5.27) 
(5.28) 



7 (x±) 

l H+H~ 
J H+H~ 



0, 



A(m 2 -m 2 H± ) 2 



+4 (m£ - m 2 H ± ) A H±k + 3 A H±k A H±l 



+D_ k D_ t 4 m\ m x ± m x ± (m x - m 2 H ± ) 
+C* k C* 2 m x m x ± [2 (m 2 x - m 2 H± ) A H±k + 3 A 2 H±k 



+C H k C» 



where A H ± k 
and ( p9| ). 



m 2 H± 



3m x ± m x ± A H±k A H±l 



2 (m x - m^±) A H±fc A#±, 



(5.29) 



m 2 x — w^±- The couplings C± k and Z)|f A are given in eqs. 



(5.30) 

(ED 



Higgs (/i, i?)-chargino (x^) interference term: 



a 



H+H- 
} H+H~ 



(5.31) 



i 7T 

fe=l 

<! c 



c 



H+H~r C^ XT 



y 

•=h,H \ X 



l k m x [2 (m- 



m^-±) + 3 A H ± 



(5.32) 



Z-chargino (Xfc) interference term: 



L H+H~ 



0: 



k=l 



c 



H+H~Z (jXXZ 



4 m 2 . — m 2 z + iTz mz 



D 



H 



(m 2 



m 



(5.33) 
(5.34) 
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4- XX -»• W+W- 

This process involves the s-channel CP-even Higgs boson (h and H) and Z exchange, and 
the t- and n-channel chargino (x^, A; = 1,2) exchange 



_ ~(h,H) ~(Z) ~( x ±) ~(fc,ff- x ±) ~(Z- X ±) 
a ww — a ww ' a WW ' a WW ' a WW ' a WW ' 

L w - + + + ^^-x±) + T(s-x±). 
www — u ww ' ww ' ww ' ww ' ww > 



(5.35) 
(5.36) 



• CP-even Higgs-boson (h, H) exchange: 



a 



(h,H) _ 
WW — u > 

r(h,H) _ 3 
W - 32 ^ 



E 



, s — + i T r m r 

r=h,H 



3 m, 



4m w ™x 



m 



w 



X . 



(5.37) 
(5.38) 



• Z-boson exchange: 



l ww 



0, 



1 



b {z) - 
°ww - 47r 



qWWZ (JW. z 



s — m z + iTz mz 
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chargino (x^) exchange: 
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+3m ± m ± (3m^ - 4m 2 ff m 2 + 4m*) A^ fc A 



16 ml m w (ml - m 2 w f A Wk A W i 
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H7 



\m w m 4 x + 4m x ) A Wk 



2 ~A~ 2m x' 

+4 m 2 x m x ± m Y ± (11 m%j - 18 m 2 ^ m\ + 16 m x ) A 2 ^ A Wi 
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where A^ 
and (p4T) . 

• Higgs (h, H )-chargino (x&) interference term: 



The couplings C^ k and D^. are given in eqs. ( 4.83; ) 
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Z-chargino (x^) interference term: 
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5. XX 



This process involves the s-channel CP-even Higgs boson (/i and .ff) exchange and the t- 
and u-channel neutralino (x® ,i = 1, . . . , 4) exchange 
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r r(h,H) . r( x °) , r{h,H- x °) 
bzz = b zz +V£ Z > + b zz * >; 

CP-even Higgs-boson (h, H) exchange: 



a 



(h,H) 
ZZ 
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64^ 
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(jZZr (jX\r 



s — mi + i T r m r 
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4. 



rrir 



neutralino (x®) exchange: 
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(5.49) 
(5.50) 



7 (X°) 
ZZ 



r(x°) 

'ZZ 



1 \^ n z*rtZ 
— Wj 



mr 



m|) 



»J=1 



16 vr ^ miAi-Ai, 

ij=l Z Z] 



xCf*Cf 3 [D§> A 2 Zl + D(? A Zi A Zj + £><f A| ?; A^ + D\f A Zi A Zj 
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where Azi 
and flPfi|) . 

• Higgs (/i, )-neutralino ) interference term: 



p. The couplings C*^ and D± { are given in eqs. Q4.95| ) 
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-m z (5 m x + 9 m o ) 



(5.58) 



- 33 - 



6- XX -»• ff 

This process involves the s-channel Higgs boson (h, H and A) and Z boson exchange and 
the t- and u-channel sfermion (f a ) exchange 
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• CP-odd Higgs-boson (^4) exchange: 
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Z-boson exchange: 



(5.59) 
(5.60) 



(5.61) 
(5.62) 



(5.63) 
(5.64) 
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where Ar 

Ja 



m 2 f — m 2 — m ~ . The index a counts sfermions so that a = 1 , . . . , 6 for squarks 



/' "7a' 

and charged sleptons, and a = 1, 2, 3 for sneutrinos. The couplings C± and -D^. are given 
in eqs. (gTo|) and ( P~l~06| ). 



Higgs (/i, i?)-sfermion (/ a ) interference term: 
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where P4 = 4 mr, — + i Ta rriA ■ 



Z-sfermion (f a ) interference term: 
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where = 4 m 2 — m 2 z + iTz mz- 



6. Summary 

The neutralino is undoubtedly the most popular candidate for a WIMP dark matter in 
the Universe. Supersymmetry remains arguably the most promising extension of the Stan- 
dard Model. The next several years will witness extensive searches for supersymmetry in 
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colliders as well as for WIMPs in underground detectors. Measurements of the cosmologi- 
cal parameters, and in particular of the relic abundance of the dark matter, have already 
reached the accuracy of a few per cent and more progress is expected. 

In light of this, theoretical computations of the neutralino relic abundance need to 
be now performed with at least the same, if not better, level of precision, if one wants 
to reliably compare theoretical predictions with observations. Motivated by this goal, we 
have derived a full set of exact, analytic expressions for the neutralino pair-annihilation 
cross sections into all tree-level two-body final states in the framework of the MSSM. 

Acknowledgments 

T.N. is grateful to Lancaster University for kind hospitality extended during his visit. The 
authors appreciate Dr. A. Birkedal-Hansen to point out several minor typographical errors 
in the first version of this paper. 



-37- 



A. The MSSM Lagrangian 



In this Appendix, the MSSM Lagrangian is given explicitly in the mass eigenstates, which 
makes it easier to read out relevant Feynman rules. 

Z-boson-Fermion-Fermion 
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and g is the usual gauge coupling of SU(2)l whereas Qf and T 3 j were defined below 
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where the matrices N, U, V were defined via eqs. ( ^ ) and (p.5|). 
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We neglect CP violation hence Cy 3 = 0. 
Gauge-Higgs-Higgs 
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Higgs-Fermion-Fermion 
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where n = 1, 2, 3 is an index for generation, so e ra = e, fj,, r, efc. 
Higgs- Gauge- Gauge 
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Note that we assume no CP violating phases in the neutralino mass matrix and we use the 
convention that TVjj is real. 
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Gauge- Gauge- Gauge 
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where 

a2 l = (r?J )L i^ L + rif )R £^ R ) , (A.49) 
Affl B = (r^L^ /n£ + r (/V rX?/ " fl ) • ( A - 5 °) 

Note that in eqs. ( |A.49 ) and ( |A.50| ), / = u, d, e, v represents the type of the fermion, 
while / in the main text denotes each fermion f = u, c, t, etc. The slepton and squark 
mass eigenstates f a (p a with a = 1,2,3 and e a ,u a and d a with a = 1,...,6) are related to 
the sfermion gauge eigenstates f n L and f n R (n = 1,2,3) via 

fnL = ^2faT a { ^ L , (A.51) 
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where Vf = (T(f)L^(f)R) denotes a 6 x 6 matrix which diagonalizes the sfermions mass 
matrix given by eqs. (|3~7|) - (^.10|) : Vf M~ VI = diag(m| , • • • , m~ ). Note that in the 

/ J1J6 

case of squarks and charged sleptons the mixing matrices rV u \ l,r, l,r an d T(e) l,r have 
dimension 6x3, while for sneutrinos T( v )l is a mixing matrix of order 3x3. (We neglect 
here the CKM matrix for simplicity.) Finally £XifnL,R anc j r x ( -fnL,R are gj ven as follows: 
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B. Auxiliary Functions 

Here we give expressions for the auxiliary functions used in the text. First, we define 

D{s,x,y 1 ,y 2 ) = x H , (B.l) 
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(Ti,yi)(s,x,yi,y 2 , zi, z 2 ) (i = 0, ■■-,4), then 



J r (s,x,y 1 ,y 2 ,z) = — In 



t+(s,x,yi,y 2 ) - z 



t-(s,x,yi,y 2 ) - z 



(B.3) 



and 



% = — - — [F(s,x,y 1 ,y 2 ,zi) - J r (s,x,y 1 ,y 2 , z 2 )], (B.4) 
z\ — z 2 

Ti = [z\J r {s,x,yi,y 2 ,zi) - z 2 J 7 (s,x,y 1 ,y 2 ,z 2 )], (B.5) 

z\ — z 2 

T 2 = l-\ [z%J 7 (s,x,y 1 ,y 2 ,z 1 ) - z%J 7 (s,x,y 1 ,y 2 ,z 2 )], (B.6) 

z\ - z 2 

T 3 = D + zi + z 2 -\ \z\ J r (s,x,y 1 ,y 2 ,z 1 ) - z\ J r (s,x,y 1 ,y 2 , z 2 )], (B.7) 

z\ — z 2 

% = D 2 + D ( Zl + z 2 ) + {zj + z\ + Zl z 2 ) + ^F 2 

H — [ziF{s,x,yi,y 2 ,zi) - z 2 F{s, x, y 1: y 2 , z 2 )], (B.8) 

zi - z 2 

^ = ; ^[-^(s,^, 2/i,2/2,^i) +F(s,x,y 1 ,y 2 ,z 2 )], (B.9) 

Zl + Z 2 — I L) 

^ = 1 — " — ^[ 2 ( Zl ~ F>)J r (s,x,y 1 ,y 2 ,zi) -2(z 2 - D) J r (s,x,y 1 ,y 2 , z 2 )], 

Z\ + Z 2 — ID 

(B.10) 

^ 2 = 1 H 1 ^tf^ 1 ( Zl ~ 2D ) :F ( s , x ,yi,y2,zi) 

z\ + z 2 -2D 

+z 2 (z 2 - 2 £>) ^(a, x, yi, y 2 , z 2 )], (B.ll) 
y 3 = 2 (zi - z 2 ) + , 1 n [21 (2 z? - 6 zi D + 4 D 2 )F(s, x, y u y 2 , z ± ) 

Z\ + Z 2 — IL> 

- z 2 (2 zl - 6 z 2 D + 4 D 2 )F(s, x, yi, y 2 , z 2 )], (B.12) 
3^4 = -# 2 - £> (zi + z 2 ) + (z 2 + z 2 - Zl z 2 ) + ^F 2 

+ - 1 on (2? - 4 zi D + 4 D 2 )F(s, x, yi , y 2 , zi) 
Zl + z 2 — Z 17 

+ z 2 (z 2 - 4 z 2 D + 4 D 2 )J 7 (s, x, yi, y 2 , z 2 )]. (B.13) 
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